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ABSTRACT 

We use a pressure-based model for splitting cold hydrogen into its atomic (HI) and molecular (H2) com- 
ponents to tackle the co-evolution of HI, H2, and star formation rates (SFR) in ~3 • 10^ simulated galaxies 
in the Millennium Simulation. The main prediction is that galaxies contained similar amounts of HI at red- 
shift 2: « 1 — 5 than today, but substantially more H2, in quantitative agreement with the strong molecular 
line emission already detected in a few high-redshift galaxies and approximately consistent with inferences 
from studies of the damped Lyman-a absorbers seen in the spectra of quasars. The cosmic H2/HI-ratio is 
predicted to evolve monotonically as fiHa/^^Hi oc (1 -I- zY'^. This decline of the H2/HI-ratio as a function 
of cosmic time is driven by the growth of galactic disks and the progressive reduction of the mean cold gas 
pressure. Finally, a comparison between the evolutions of HI, H2, and SFRs reveals two distinct cosmic 
epochs of star formation: an early epoch (z > 3), driven by the evolution of Qiii^ii^{z), and a late epoch 
(z < 3), driven by the evolution of flii^{z)/fliii{z). 

Subject headings: galaxies: high-redshift — galaxies: evolution — ISM: atoms — ISM: molecules — cosmol- 
ogy: theory 



1. INTRODUCTION AND KEY IDEA 



JL Neutral hydrogen is the fuel for the formation of stars, 
^^he cosmic star formation rate (SFR) density as inferred 
"*Tfrom ultraviolet, far-infrared, and submillimeter observa- 
j^ipns incr eases by an ord er of magnitude from redshift z = 
I— tti z = 2 ()Hopking||2007D . Hence, neutral hydrogen in early 
.g alaxies was either more abundant or transformed into stars 
''^riore efficiently than today. 

A useful quantity in this context is the star formation 

~^ efficiency (SEE) of a galaxy, defined as the SFR divided 

the gas mass. The weak cosmic evolution of the den- 
(-~§ity of neutral atomic hydrogen (HI), deri ved from Lyman- 
. alpha absorptio n against distant quasars ()Lah et al.l l2007t 
">!|P6ntzen fc Pett ini 2009), indicates a strongly increased SEE 
^l^t high z. But recent detections of str ong molecular line 
Sipmission in ordinary galaxies at z = 1.5 ijDaddi et al.ll2008l ) 
T^4gg6st that the SFEs of these galaxies are similar to those 
K^seen today. The seeming contradiction between these two 
'^^^onclusions arises from the conceptual confusion of SFEs 
r^nferred from galactic HI with those inferred from H2. In 
i3act, it is crucial to distinguish between the two quantities 
■ ■ SFEhi = SER/A/hi and SFEh^ = SFR/Mh^- In principle, 
there is no contradiction between the detected strong cosmic 
evolution of SFEhi and the weak evolution of SFEhj - these 
empirical findings could simply imply that the H2/HI-mass 

ratios ii^f^f^ of galaxies increase substantially with z. 

In this letter, we show that there is indeed strong theoret- 
ical support for such an increase of R^^^i^^ with z in regular 
galaxies. This evolution is driven by the approximate scal- 
ing of galaxy sizes as (1 + z)~^ predicted by dark matter 
theorv (jGunn &: Gottlll97§) and confirnied by observations 
in the Ultra Deep Field (|Bouwens et"al]|2004j) . Hence, the 
cold gas disks at high redshift must, on average, be denser 
than today. Combining this prediction with the relation 
bet ween gas pressure and H2/ HI-ratios in nearby galaxies 
fe.g. lBlitz fc Rosolowskvl[20060 . leads to the conclusion that 

^moi^^ must increase dramatically with z. Our quantita- 
tive predictions of this evolution rely on a recently presented 
semi-analytic numerical simulation of HI and H2 in ^3 • 10^ 



simulated galaxies (lObreschkow et all I2009D . based on the 
Millennium Simulation (Springel et al. 2005). 

Section [2] overviews our simulation method and the model 
for the H2/HI-ratio in galaxies. In Section [3l we present and 
interpret the predicted evolution of galactic HI and H2 and 
their relation to star formation. Section |4] compares these 
predictions to empirical data, and Section [5] summarizes our 
key conclusions. 

2. SIMULATING HI AND H2 IN GALAXIES 
2.1. Physical model for galactic H2/HI-ratios 
In virtually all re gular galaxies in the local Universe, 



whether spir als (e.g. iLerov et al.l l2008l ) or ellipticals (e.g. 
lYoungj [2002! ) ■ the cold gas resi des in a flat disk. Some 
observations of CO at 2; « 2 (jTacconi et al.l I2006D sug- 
gest that even at high redshift most cold gas lies in disks. 
Based on this evidence, we have recently introduced a 
model for the dis tributions of HI and H2 in regular galaxies 
(jObreschkow et a l. 2009) , assuming that all cold gas resides 
in a flat symmetric disk with an exponential surface density 
proflle and that the l ocal H2/HI-ratio is dictated by the kine- 
matic gas pressure ([Blitz fc Rosolowskvl 120061 : iLerov et"aLl 
l2008f l. Within these assumptions, we could show that the 
H2/HI-mass ratio R^^^i^^ of an entire galaxy is given by 
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where R'^^i represents the H2/HI-ratio at the galaxy center, 
can be approximated as 



R'^,, = [ll.3m4kg-2r-4Mg(Mg-|-0.4Mf ^i^)] 



0.8 



(2) 



where r^isk is the exponential scale radius of the disk, Mg 
is the total cold gas mass, and AI^^^^ is the stellar mass 
in the disk. Eqs. (|ll2p constitute a physical model to esti- 
mate Rf^li^^ in regular galaxies based on AI^^^^, Mg, and 



galaxy 
mol 



Tdisk ■ In order to predict the cosmic evolution of R 
therefore require a model for the co-evolution of M^^^^, Mg, 
and Tdisk in galaxies. To this end, we adopted the virtual 
galaxy catalog of the Millennium Simulation described in 
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Section [2T2l The limitations of the model of Eqs. (|ll2p and 
their impact on the pred icted H2/HI-ratios are discussed in 
(|Obreschkow et al.|[2009l) . 

2.2. HI and H2 in the Millennium Simulation 

The Millennium Simulation fS pringel et al.ll2005[ ) is an TV- 
body simulation within the ACDM cosmology of ~10"'^° grav- 
itationally interacting particles in a periodic box of comoving 
volume (500 Mpc)^, where Hq = 100 /ikms^^ Mpc^^ 
and h = 0.73. The evolving large-scale structure generated 
by this simulation served as the skeleton for the simulation of 
~3 • 10^ galaxies at th e halo centers. In the "sem i-analytic" 
approach adopted by iDe Lucia fc Blaizo^ (|2007l ). galaxies 
were considered as simplistic objects with a few global prop- 
erties that are evolved stepwisely using a list of physical pre- 
scriptions. For example, the total amount of cold hydrogen 
(HI+H2) in a galaxy is defined by the history of the net ac- 
cretion, which in the model consists of (i) the infall of gas 
from the hot halo, (ii) the loss of gas by star formation, and 
(iii) outflows driven by supernovae and active galactic nuclei. 
Star formation in each galaxy is tackled using a law, where 
all cold gas above a critical surface density is transformed 
into stars on a timescale p roportional to the dynamical time 
of th e disk (for details s e elCrqton et~al. 20061. 

In I Obreschkow et al] (|2009D . we applied the model of 
Section 12.11 to the simula ted galaxies in the catalog of 
IDe Lucia fc BlaizotI (|2007( ) ("DeLucia-catalog"), to split 
their cold hydrogen masses into HI and H2. Our simulation 
successfully reproduced many local observations of HI and 
H2, such as mass functions (MFs), mass-diameter relations, 
and mass-velocity relations. Yet, the high-redshift predic- 
tions are inevitably limited by the semi-analytic recipes of 
the DeLucia-catalog. The most uncertain recipes are those 
related to mergers (e.g. feedback of black hole coalescence 
and starbursts), but they have a minor effect on the cosmic 
space densities of HI and H2 , since most cold gas in the sim- 
ulation is found in regular disk galaxies^ with at most minor 
merger histories. However, inaccurate prescriptions for iso- 
lated galaxies could significantly affect the space densities of 
HI and H2, and it may well become necessary to refine our 
simulation as improved semi-analytic methods come on line. 

3. RESULTS 
3.1. Predicted evolution of HI and H2 

Fig. [1] shows the predicted evolution of the HI-MF and 
H2-MF, i.e. the comoving space densities of sources per log- 
arithmic mass interval. The predictions at z = roughly 
agree with available observational data, but the obvious dif- 
ferences, such as the spurious bumps around Mm ~ 10*'^ 
and Mho ~ 10^ (a mass re solution limit), have been dis- 
cussed in lObreschkow et "all (|2009( ). 

The predicted Hl-masses remain roughly constant from 
z = to 2: — 2, while H2-masses increase dramatically. 
These different evolutions are also reflected in the comov- 
ing space densities ilni = Pm/Pc and flu^ = PH2/PC, where 
Pc{z) = 3iJ^(z)/(87rG) is the critical density for closure. 
Here, ftm and r^Ha only account for gas in galaxies, ex- 
cludi ng unbound HI between th e first galaxies (iBecker et al.l 
|2001[ ) or possible H2 in haloes (jPfenniger fc Combeslll994[) . 
The simulated functions flmiz) and Q}i^(z) are shown in 
Figs. [2^, b, while Fig.[2t represents their ratio i?™™"^(2) = 

^ By contrast, a significant fraction of the stars at z = is in massive 
elliptical galaxies with violent merger histories, but even those galaxies 
formed most stars in their spiral progenitors. 



il}i^{z)/Q,-iiiiz), which is closely described by the power-law 

i?~(z) « 0.3 • (1 + 2)1-^ (3) 

The simulation yields R^^^^^'^iO) « 0.3 and finds the 
crossover, i?™™"^(z) = 1, at z « 1.4. Our model pre- 
dicts that Eq. ^ extends to epochs, where the first galaxies 
formed, but this prediction is likely to breakdown at the 
highest redshifts, whe re the formation of H 2 was inhibited 
by the lack of metals (jAbel fc Haimanll2000l ). 
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Fig. 1. — MFs of HI and H2. Lines show the simulation results 
at z = (solid), z = 2 (dashed), z = 5 (dash-dotted), z = 
10 (dotted). Squa re dots represent the empirical data and 1-cr 
scatter at z = (|Zwaan et al.l l2005al : lObreschkow fc Rawlingj 
12009131 ). and the open circle represents our den sity estimate at 
z = 1.5 (Section [H) based on lDaddi et~all l|2008l '). 

Physically, the strong evolution of H2/HI is essentially 
driven by the size-evolution of galaxies and their haloes. 
The Millennium Simulation assumes that the virial radius 
fvir of a spherical halo always encloses a mass with an av- 
erage density 200-ti mes above the critical density oc H^ 
(jCroton et al.ll200^ . Hence, for a fixed halo mass, r^ir oc 
H~'^/^ . In a flat Universe this implies 

rvirOC [an(l + z)'+f^A]"'^', (4) 

which asymptotically tends to ryjr oc {\ + z)~^ f or high z. By 
virtue of the theory of iFall fc Efstathioul (|1980[ ). this cosmic 
scaling of rvir results in a similar scaling of the disk radius, 
i.e. Tdisk oc (1 -|- z)~^ , consistent with observations in the 
Uhra Deep Field (jBouwens et al.ll2004.1 . 

For the gas-dominated galaxies in the early Universe, 
Eq. (HI) reduces to Rl^^^ oc r^^^^^M^-^. Yet, the cold gas 
masses Mg of individual galaxies in the simulation evolve 
weakly with cosmic time, due to a self-regulated equilib- 
rium between the net inflow of gas and star formation. In 
fact, most of the evolution of Jlni+Ha in the redshift range 
z « 3 — 10 is due to the build-up of new galaxies. There- 
fore, Rl^^^ oc r~l^ oc (1 -h zf-^. At redshifts z « 1 - 10, 
-^moi typically takes values between 10 and 10'*, such that 
Eq. ^ can be approximated as R^^i^^ -^moi°^- Hence, 
^moi^^ OC (1 -I- z)^'^, which explains the scahng of Eq. 

The cosmic evolution of l^Ha shown in Fig. [2] can be di- 
vided in two epochs: The early epoch (z > 3), where fln^ 
increases with cosmic time, and the late epoch (z < 3), where 
S1h2 decreases with time. In the early epoch, Rf^li^^ > 1 
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Fig. 2. — Cosmic evolution of the fractional space densities of 
HI and II2. Solid lines represent the simulated evolution of f^Hi 
(a), (b), and r^Ha/^^Hi (c). The dashed line in panel (c) is the 
power- law fit for Q,ii^/^m given in Eq. (|3J. The points represent 
the observations described in Section U) The dash-dotted lines 
in panels (a) and (b) represent the evolution of f2Hi and r^Hj, 
if the correction factor ^ in Obrcschko w et alj (|200 9;') is set to 
(" = 1, i.e. the total cold gas mass of the DeLucia-catalog is not 
corrected. We originally introduced this correction to fit the sum 
of the local space densities of HI and H2 (shown as diamonds). 

implies ^Yi2 ~ ^hi+H2 i hence the growth of f^Ha re- 
flects the general increase of fini-i-Ha due to the intense as- 
sembly of new galaxies. In the late epoch, R^^^^ < 1 im- 
plies that r^Ha ~ -Rmo™"^^HH-H2- At this epoch the forma- 
tion of the massive galaxies in the simulation is completed, 
i.e. VLm+U2{z) ~ const and f^Ha oc i?™^?"'". Thus the de- 
crease of r^Ha ill this late epoch is driven by cosmic decline 
in i?™^™"^ or, physically, by the cosmic evolution of pressure. 

3.2. Link between HI, H2, and star formation 

To discuss the global cosmic evolution of the efficiencies 
SFEhi and SFEhs (Section [T]), we shall define 

(SFEhi) = Psfr/phi , (SFEH2) = Psfr/ph2, (5) 

where pm oc 0.m, PH2 Vl]^^ , and psFR denote the comoving 
space densities of HI, H2, and SFR. 

In the semi-analytic recipes of the DeLucia-catalog, SFRs 
are estimated from the gas density and the dynamical time 
scale of the disk (Section|2.2l) . This Schmidt-Kennicutt law 
(|Schmidti Il959t iKennicutt 19983 for star formation makes 
similar predictions to mod els based on cold gas pressure(e.g. 
iBlitz fc Rosolowskvl [20061 ) . and therefore the SFRs in the 
DeLucia-catalog are, by default, approximately consistent 
with our model to split cold hydrogen into HI and H2. The 
evolutions of (SFEhi) and (SFEH2) predicted by the simu- 
lation again reflect the marked difference between HI and 
H2. They are approximated (~20% relative error) by the 
power-laws, 

(SFEHi)/[G2/r-i] = 0.23 (1 + zf-\ (6) 
(SFEH2>/[Gyr-i] = 0.75 (1 + (7) 
out to z « 8. 
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Fig. 3. — A simplistic model for the cosmic history of star formation. 



Due to the low power in Eq. ([7]) psfk{z) is approximately 
proportional to n]i2{z). We can therefore apply the two 
cosmic epochs of r^Hs {z) introduced in Section 13.11 to the 
history of star formation (see Fig. [3]): In the early epoch 
(z ^ 3), psFR increases with cosmic time, proportionally 
to riHi+H2- This increase traces the dramatic assembly of 
new galaxies. In the late epoch {z < 3), psFR decreases 
roughly proportionally to r2H2/^Hi- This epoch is driven 
by the cosmic evolution of pressure (or density) in galactic 
disks. This interpretation of the history of star formation 
does not, in fact, conflict with the picture that star forma- 
tion is ultimately defined by the accreted cold gas mass (see 
Section [2.2p and a Schmidt-Kennicutt law for transforming 
this gas into stars. Our H2/HI-based interpretation simply 
adds another layer to the causal chain, by suggesting that 
cold gas mass and density ultimately dictate the amount of 
molecular material available for star formation. 

The simulation also includes star formation via merger- 
driven starbursts, associated with the creation of the stellar 
spheroids of early- type spiral or elliptical galaxies. However, 
the cosmic star formation density caused by mergers only ac- 
counts for about 1% of psFR in the semi-analytic simulation 
of the DeLucia-catalog. (This fraction should not be con- 
fused with the fraction of gas- rich or "wet" mergers, since 
only some of the cold gas involved in these mergers is effi- 
ciently turned into stars in the form of a starburst.) 

4. COMPARISON WITH OBSERVATIONS 

The DeLucia-catalog and our post-processing to assign 
HI and H2, rely on established data of the local Universe. 
Our simulated HI- and H2-properties at z = are consis- 
tent with all available observatio ns, i.e. MFs (see Fig. [1] ), 
disk sizes, and velocity profiles (jObreschkow et aP l2009f i . 
In particular, the simulated values r2Hi(0) = 3.5 ■ 10~^ 
and Vl-n^ifS) = 1.2 • 10~* are consistent with the val- 
ues (diamonds in Fig. [2]) derived fro m the MFs observed 
in HI- and CO-emission at z w (jZ waan "eFaH [2005at 
lObreschkow fc RawHngsl l2009b| ) . At z > 0, the currently 
available data are sparse, especially in emission. 

The only measurement of flm in emission at intermediate 
redshift is based on the stacking of 121 galaxies dX z — 0.24 
(jLah et a l. 2007, tr iangle in Fig. ^ The detection is specu- 
lative (see Fig. 7 in lLah et al.ll2007t ), but roughly consistent 
with our simulation. All other measurements of J^hi at z > 
rely on absorption detections of dam ped Lyman-g sy stems 
(DLAs) . Respectiv e data points fro m Rao e t all (|2006 ^ (cir- 
cles in Fig. [2]) and iProchaska et a l. (2005) (filled squares) 
are, taken together, inconsistent with the predicted values 
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of r^Hi- By contrast, IZwaan et al.l (|2005H ) demonstrated 
that the population of Hl-galaxies in the local Universe can 
fully explain the column density distributions of DLAs out 
to z = 1.5, consistent with the nearly absent evolution of 
r^Hi from z = 0tO2; = 1.5 predicted by our simulation. At 
present it is therefore difficult to judge, whether the simu- 
lation is inconsistent with empirical data at these low red- 
shifts covering 2/3 of the age of the Universe. At higher 
redshifts, however, the measurements of JIhi seem not rec- 
oncilable with the simulated result, and even accounting for 
gravitational lensing by the DLAs only corrects the empir- 
ical values of I^hi by about 30% (open squares in Fig. [21 
iProchaska et al.l l2005l ) . The simulated values of fini are 
likely to underestimate the real values by about a factor 
2 - a plausible offset given the long list of simplifying ap- 
proximations required from the iV-body Millennium Simu- 
lation to our final post-processing of hydrogen in galaxies. 
Much progress could be expected from treating Hl-masses 
and H2-masses as separate quantities directly in the semi- 
analytic galaxy simulation. This would allow, for example, 
to refine the feed back-mechanisms fo r suppression of gas in- 
fall (explained in lCroton et al.|[2006( ). such that HI can still 
be accreted, while the formation of H2 and stars is inhibited. 
Such a semi-analytic setting would also allow the implemen- 
tation of a recipe for the large-scale dissociatio n of molecu- 
lar g as by the radiation of newly formed stars (lAllen et al.l 
Il986l ). Both examples would effectively increase the amount 
of HI in high-redshift galaxies. 

The most representative high-redshift observations of 
molecular gas to-date rely on two normal galaxies [BzK- 
4171 and BzK-2lQ m) ai z k. 1.5, r eUably detected in 
C0(2-l) emission bv lDaddi et all (|2008f) . Unlike other CO- 
sources at similar or higher z, these objects are ordinary 
massive galaxies with FIR-luminosities of ipiR ~ 10^^ Lq, 
selected only due to the availability of precise spectroscopic 
redshifts. From these two detections, we estimated the 
H2-space density (empty circle in Fig. [T]) as follows: The 
mass interval spans between the masses « 2 ■ lO^^M© 
and ~ 3 ■ 10^°, respectively obtained for EzK-Alll 

and BzK-2\QQQ by applying th e CO-to-H2 conver sion of 
a = 1 M0(K km s'lpc-^)-! (jPaddi et al.l[200l . The 
space density of these CO-sources was approximated as the 
space density of FIR-sources at Lfir ~ 10^^ Lq, based 
on the fact that all (both) targeted galaxies with Lfir ~ 
10^^ L0 revealed similar CO-luminosities Leo- We esti- 
mate their space density to be 1 — 2 • 10~^ Mpc~'^ per unit 



of log(L FiR.) by extrapolatin g the FIR-luminosity functions 
(LFs) of iHuvnh erall (|2007l ). Since Lpm oc Leo « Mh^, 
we find roughly the same space density per unit log(MH2)- 
These result is consistent with the simulated H2-MF a± z — 2 
(Fig. [J). 

C onsidering H2-absorption studies, ICurran et al.l (|2004f ) 
and lNoterdaeme et al.l (|2008f ) have determined H2/HI-ratios 
in DLAs that showed H2-absorption. They found H2/HI- 
ratios of ~ 10~^ to ^ 10^^ at z « 2 — 3, clearly much 
smaller than our prediction for Vt-n.^/Vlni. We argue that 
measurements of H2/HI in DLAs do not trace riH2/^Hi 
since DLAs are by definition HLselected objects and H2 has 
a much smaller space coverage than HL In fact, H2-disks 
in gala xies are much smaller than Hl -disks, especially at 
high z (jObreschkow fc Raw lings"2009a1 , and even inside the 
H2-disks the coverage of H2 is small compared to HI (e.g. 
lFerrierell200ll) . A more detailed explanation of why H2- 
searches in DLAs ar e expe cted to be difficult was given by 
IZwaan &: Prochaskal (|2006l ) based on the analysis of CO- 
emission maps of local galaxies. 

5. CONCLUSIONS 

In this letter, we have predicted the cosmic evolution of 
HI- and H2-masses in ~ 3 • 10^ simulated galaxies based on 
the Millennium Simulation. The predicted cosmic decline in 
the H2 /Hl-ratio is consistent with the weak cosmic evolution 
of r^Hi inferred from DLA-studies and recent observations 
revealin g a significantly e nhanced space density of H2 at 
z = 1.5 (|Daddi et al.ll200l . 

Perhaps the most important conclusion is that HI- and H2- 
masses evolve very differently with cosmic time and there- 
fore cannot be used as proportional tracers of one another, 
especially not for the purpose of high-redshift predictions. 
There is no contradiction between the large H2-masses de- 
tected at high z, which imply values of SFEna similar to 
those in the local Universe, and the weak evolution of HI, 
implying massively increased values of SFEhi at high z. 
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